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Exotic dark matter together with dark energy or cosmological constant 
seem to dominate in the Universe. An even higher density of such mat- 
ter seems to be gravitationally trapped in our Galaxy. The nature of dark 
matter can be unveiled only, if it is detected in the laboratory. Thus the 
accomplishment of this task is central to physics and cosmology. Current 
fashionable supersymmetric models provide a natural dark matter candi- 
date, which is the lightest supersymmetric particle (LSP). Such models 
combined with fairly well understood physics like the quark substructure 
of the nucleon can yield the effective transition operator at the nucleon 
level. We know, however, that the LSP is much heavier than the proton, 
while its average energy is in the keV region. Thus the most likely possibil- 
ity for its direct detection is via its elastic scattering with a nuclear target 
by observing the recoiling nucleus. In order to evaluate the event rate one 
needs the nuclear structure (form factor and/or spin response function) 
for the special nuclear targets of experimental interest. Since the expected 
rates for neutralino-nucleus scattering are expected to be small, one should 
exploit all the characteristic signatures of this reaction. Such are: (i) In the 
standard recoil measurements the modulation of the event rate due to the 
Earth's motion, (ii) In directional recoil experiments the correlation of the 
event rate with the sun's motion. One now has both modulation, which is 
much larger and depends not only on time, but on the direction of obser- 
vation as well, and a large forward-backward asymmetry, (iii) In non recoil 
experiments gamma rays following the decay of excited states populated 
during the Nucleus-LSP collision. Branching ratios of about 6 percent are 
possible, (iv) novel experiments in which one observes the electrons prod 
used during the collision of the LSP with the nucleus. Branching ratios of 
about 10 per cent are possible. 
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1. Introduction 

One of the greatest discoveries of the last couple of decays is the realiza- 
tion that the matter, which can shine, constitutes a very small fraction of 
the universe. The universe is composed primarily of dark matter and dark 
energy, which are incapable of emitting light. If one defines fij = ^ for the 
species i, where p c the critical density (see sec. 2), the recent cosmological 
observations [1], when combined, lead to: 

tt b = 0.05, £Icdm = 0.30, Q A = 0.65 

for ordinary matter, dark matter and dark energy respectively. Their sum 
is [1] 17 = 1.11 ± 0.07, i.e. unity within the experimental errors , which 
means that the universe is flat. 

Additional indirect information comes from the rotational curves [2]. 
The rotational velocity of an object increases so long is surrounded but 
matter. Once outside matter the velocity of rotation drops as the square 
root of the distance. Such observations are not possible in our own galaxy. 
The observations of other galaxies, similar to our own, indicate that the 
rotational velocities of objects outside the luminous matter do not drop 
(see Fig. 1). So there must be a halo of dark matter out there. 

As we have already mentioned, however, it is only the direct detec- 
tion of dark matter in the laboratory, which will unravel the nature of the 
constituents of dark matter. In fact one such experiment, the DAMA, has 
claimed the observation of such signals, which with better statistics has 
subsequently been interpreted as modulation signals [3]. These data, how- 
ever, if they are due to the coherent process, are not consistent with other 
recent experiments, see e.g. EDELWEISS and CDMS [4]. 
Supersymmetry naturally provides candidates for the dark matter con- 
stituents. In the most favored scenario of supersymmetry the LSP can be 
simply described as a Majorana fermion (LSP or neutralino), a linear com- 
bination of the neutral components of the gauginos and higgsinos [6]- [9]. 
We are not going to address issues related to SUSY in this paper. Most 
SUSY models predict nucleon cross sections much smaller the the present 
experimental limit as < 10~ 5 pb for the coherent process. As we shall see 
below the constraint on the spin cross-sections is less stringent. 

Since the neutralino is expected to be non relativistic with average ki- 
netic energy < T >« AOKeV (m x /100GeV) , it can be directly detected 
mainly via the recoiling of a nucleus (A,Z) in elastic scattering. 

In some rare instances the low lying excited states may also be populated 
[11]. In this case one may observe the emitted 7 rays. Finally it has also 
recently been suggested [23] that one may be able to detect the neutralino 
by observing the ionization electrons produced during its collision with the 
nucleus. 
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M33 rotation curve 



Figure 1. The rotational velocity distance. Clearly it does not fall outside the luminous 
matter. This is suggestive of the presence of a dark matter halo outside the M33 galaxy. 



In every case to extract from the data information about SUSY from the 
relevant nucleon cross section, one must know the relevant nuclear matrix 
elements [12]- [13]. The static spin matrix elements used in the present work 
have been obtained from the literature [11]. 

Anyway since in all types of experiments the obtained rates are very low, 
one would like to be able to exploit the modulation of the event rates due 
to the earth's revolution around the sun [14], [15]. In order to accomplish 
this one adopts a folding procedure, i.e one has to assume some velocity 
distribution [14, 15, 17, 18] for the LSP. One also would like to exploit 
other signatures expected to show up in directional experiments [19]. This 
is possible, since the sun is moving with relatively high velocity with respect 
to the center of the galaxy. 

2. Matter Density Versus Cosmological Constant (Dark Energy). 

As we have mentioned in the previous section the universe is primarily 
made of dark matter and dark energy. The nature of each one of them 
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is not presently known. We hope the future experiments to shed light on 
this. The future of the universe depends on the composition and the nature 
of these two constituents. In the present case we will briefly explore the 
possibility that the dark energy density is constant, independent of the size 
of the universe, i.e. it coincides with Einstein's cosmological constant. 

The evolution of the Universe is governed by the General Theory of 
Relativity. The most commonly used model is that of Friedman, which 
utilizes the Robertson- Walker metric 

(ds) 2 = (dtf - a 2 (t)[-^^ + r 2 {{d6) 2 + sin 2 #(#) 2 )] (1) 

, where a(t) is the scale factor. The resulting Einstein equations are: 

R^u ~ \g»uR = SitGnT^ + k 9lxv (2) 

where Gat is Newton's constant and A is the cosmological constant. The 
equation for the scale factor a(t) becomes: 

d 2 a 4vr AnG N p A 

^ = -Y GNip + 3p)a = 3— ° + 3 (3) 

where p is the mass density. Then the energy is 

1 ar 2 „ rn , A 2 k 
h = — ?Ti( — — ) — Gat — {Airpa ) + —ma = constant = ——m (4) 

•—i Ctu Ob D 



This can be equivalently be written as 

H 2 + ~ 2 = Y G NP+ - (5) 



where the quantity H is Hubble's constant defined by 

H= 1 -^ (6) 

adt w 

Hubble's constant is perhaps the most important parameter of cosmology. 
In fact it is not a constant but it changes with time. Its present day value 
is given by 

#o = (65 ± 15) km/s M' 1 (7) 

In other words # -1 = (1.50±).35) x 10 10 y, which is roughly equal to the 
age of the Universe. Astrophysicists conventionally write it as 

H = 100 h km/s M pt } , 0.5 < h < 0.8 (8) 
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Equations 3-5 coincide with those of the Newtonian theory with the fol- 
lowing two types of forces: An attractive force decreasing in absolute value 
with the scale factor (Newton) and a repulsive force increasing with the 
scale factor (Einstein): 

uiM 1 
F = —Gn — — (Newton) , F = -Ama (Einstein) (9) 

Qi o 

Historically the cosmological constant was introduced by Einstein so that 
General Relativity yields a stationary Universe, i.e. one which satisfies the 
conditions: 

Indeed for k > 0, the above equations lead to a = a c =constant provided 

*^ a ^ 14 14 

-Aa c - —GnPclc = , -Aa 2 . - —Gat pa 2 = k (11) 

These equations have a non trivial solution provided that the density p and 
the cosmological constant A are related, i.e. 

A = AttG nP (12) 

The radius of the Universe then is given by 

-tj^r < i3 > 

In our days we do not believe that the universe is static. So the two pa- 
rameters, density and cosmological constant, need not be related. Define 
now 

= — , = — , p v = — — — (" vacuum" density) (14) 
Pc Pc 8itGn 

The critical density is 

9 o , 9 q n nucleons 
Pc = 1.8 x l(T 23 /i 2 ^ = 10/i 2 ^ 15 

cm 6 m 6 

With these definitions Friedman's equation E = —k 1 ^ takes the form 



K 



a 
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(n m + n A - i)h 2 (16) 



Thus we distinguish the following special cases: 

k>0 44> fl m + J7a > 1 o Closed curved Universe (17) 
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k = 44> Q. m + f^A = 1 "v^ Open Flat Universe (18) 

k < 44> f2 m + < 1 Open Curved Universe (19) 

In other words it is the combination of matter and "vacuum" energy, which 
determines the fate of the our Universe. 

To proceed further and obtain a(t) we need an equation of state, which 
relates the pressure and the density. This of the form: 

p = wpc 2 

the one finds 

a(t) = const ^(3(1+-)) , £W = _^E {1 + 3w)p (20) 

We thus see that we have acceleration if w < — 1/3 and deceleration other- 
wise. The following cases are of interest: 

— w = 1/3. This corresponds to a universe dominated by radiation 

— w = 0. This corresponds to a universe dominated by non relativistic 
matter. 

— w = — 1. This means that the cosmological constant is dominant. In 
this case the Hubble parameter is a constant and a{t) = aoe Hot . 

— — l<w< —1/3. This corresponds to dark energy. 

We are not going to elaborate further on this issue. The observational sit- 
uation at present is summarized in Fig. 2. 



3. Rates 

The differential non directional rate can be written as 

pCO^) in 

dRundir = da(u,v)\v\ (21) 

where da(u, v) was given above, p(0) = O.SGeV/cm 3 is the LSP density in 
our vicinity, m is the detector mass and m x is the LSP mass 

The directional differential rate, in the direction e of the recoiling nu- 
cleus, is: 



dR dir = ^--^—\v\v.e 0(v.e) — da(u, v S( v - v.e) (22) 
where Q(x) is the Heaviside function and: 

da(u, v) == 2 ^ U bv)2 [&sF{u) 2 + ZspinFu («)] (23) 
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Figure 2. The allowed region in the fl m , 0.de, with m indicating all matter and DE 
the dark energy, for the various experiments indicated on the plot. 



where u the energy transfer Q in dimensionless units given by 
« = 7T > Qo = [m p Ab}- 2 = 40 A' 4 / 3 MeV 



(24) 
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TABLE 1. The static spin matrix elements for various nuclei. For 
light nuclei the calculations are from DIVARI. For 127 1 the results 
presented are from Ressel and Dean (*) and the Finish group (**). 
For 207 Pb they were obtained by the Ioannina team (+). 





19 p 


29 Si 


23 Na 127 7 * 


127 j** 


207 Pb+ 



[fio(O)] 2 


2.610 


0.207 


0.477 


3.293 


1.488 


0.305 


2.807 


0.219 


0.346 


1.220 


1.513 


0.231 


f2 (0)fii(0) 


2.707 


-0.213 


0.406 


2.008 


1.501 


-0.266 


Hth 


2.91 


-0.50 


2.22 










2.62 


-0.56 


2.22 










0.91 


0.99 


0.57 









with b is the nuclear (harmonic oscillator) size parameter. F(u) is the nu- 
clear form factor and Fu(u) is the spin response function associated with 
the isovector channel. It is the area of constructing such form factors and 
momentum distributions that Professor M. Grypeos, honored in this vol- 
ume, has made a great contribution [20] . 
The scalar cross section is given by: 



a P,x oA 



1 + fh 2Z-A 

1 _!_ fs 
1 + 7<T 

Js 



a 



N, X ' 



o( 



fj, r {p)' 



,2 A 2 



(25) 



(since the heavy quarks dominate the isovector contribution is negligible). 

a% o is the LSP-nucleon scalar cross section. The spin Cross section is 

■ x u 
given by: 



J spin 



= ( 



\ 2 spin 

' a P,x° 



Cspini Cspin — 



3(1 + 4: 

J A 



-S(u) 



S(u) « S(0) = [(^-^o(O)) 2 + 2^fi (0)f2i(0) + fti(0)) s 

J A J A 



(26) 



(27) 



Some static of interest to the planned experiments are given in table 1. The 
shown results are obtained from DIVARI [13], Ressel and Dean (*) [12], the 
Finish group (**) [16] and the Ioannina team (+) [5], [10]. 

f% f\ are * ne isoscalar and the isovector axial current couplings at 
the nucleon level obtained from the corresponding ones given by the SUSY 
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models at the quark level, f A (q), f\{q), via renormalization coefficients g\, 
9a, i-e- /° = 3a /a (?) '/A = 3a/A(<?)- These couplings and the associated 
nuclear matrix elements are normalized so that, for the proton at u = 
0, yield Cspin = 1- If the nuclear contribution comes predominantly from 
protons (Q\ = S7o = Q p ), S(u) « ftp and one can extract from the data the 
proton cross section. If the nuclear contribution comes predominantly from 
neutrons ((fio = — ^1 = Q n ) one can extract the neutron cross section. 
In many cases, however, one can have contributions from both protons 
and neutrons. The situation is then complicated, but it turns out that 
g A = 0.1 , g\ = 1.2 Thus the isoscalar amplitude is suppressed, i.e 
5(0) w Slf. Then the proton and the neutron spin cross sections are the 
same. 



4. Results 

To obtain the total rates one must fold with LSP velocity and integrate 
the above expressions over the energy transfer from Q m in determined by 
the detector energy cutoff to Q ma x determined by the maximum LSP ve- 
locity (escape velocity, put in by hand in the Maxwellian distribution), i.e. 
v esc = 2.84 vq with vq the velocity of the sun around the center of the 
galaxy (229 Km/s). 

4.1. NON DIRECTIONAL RATES 

Ignoring the motion of the Earth the total non directional rate is given by 



R = K 



Ccoh{A, Hr{A))(T^ x o + C spin (A, fl r (A))a^ Cspin (28) 



where K 



P(0) 



\J {v 2 } and 



C coh (A, fi r (A)) 



Mr (A) 
fl r (p) 



^4 t co h(A) , C sp i n (A, fl r (A)) 



fir (A) 



-1 2 



fj, r (p) 



spin 



(A) 



A 



where t is the modification of the total rate due to the folding and nuclear 
structure effects. It depends on Q m i n , i.e. the energy transfer cutoff im- 
posed by the detector and a = [fi r bvo V^]" 1 The parameters c co h{A, fj, r (A)), 
c S pin(A, fi r (A)), which give the relative merit for the coherent and the spin 
contributions in the case of a nuclear target compared to those of the pro- 
ton, are tabulated in table 2 for energy cutoff Q m in = 0, 10 keV. Via Eq. 
(28) we can extract the nucleon cross section from the data. 



10 



TABLE 2. The factors cl9 = c co/l (19, M 19 )); sl9 = c spi „(19, /i r (19)) and 
cl27 = c coh (127,Mr(127)), s!27 = c spin (127, ^(127)) for two values of Q min . 



Qmin 










m x 


(GeV) 








keV 




20 


30 


40 


50 


60 


80 


100 


200 





cl9 


2080 


2943 


3589 


4083 


4471 


5037 


5428 


6360 





sl9 


5.7 


8.0 


9.7 


10.9 


11.9 


13.4 


14.4 


16.7 





cl27 


37294 


63142 


84764 


101539 


114295 


131580 


142290 


162945 





sl27 


2.2 


3.7 


4.9 


5.8 


6.5 


7.6 


8.4 


10.4 


10 


cl9 


636 


1314 


1865 


2302 


2639 


3181 


3487 


4419 


10 


sl9 


1.7 


3.5 


4.9 


6.0 


6.9 


8.3 


9.1 


11.4 


10 


cl27 





11660 


24080 


36243 


45648 


58534 


69545 


83823 


10 


sl27 





0.6 


1.3 


1.9 


2.5 


3.3 


4.0 


5.8 



Using fif 



1.22 and Of 



2.8 for 127 I and 19 F respectively the extracted 



nucleon cross sections satisfy: 



a 



spin 
P,X° 

a P ,x° 



CcohjA, (J, r (A)) 
{A,ll r {A)) 



=>« xlO 4 (A 

1 La 



127) , « xlO 2 (A = 19) 
(30) 

It is for this reason that the limit on the spin proton cross section extracted 
from both targets is much poorer. For heavy LSP, > 100 GeV, due to the 
nuclear form factor, t sp i n (127) < t sp i n (19). This disadvantage cannot be 
overcome by the larger reduced mass. It even becomes worse, if the effect 
of the spin ME is included. For the coherent process, however, the light 
nucleus is no match (see Table 2). 

If the effects of the motion of the Earth around the sun are included, the 
total non directional rate is given by 



R = K 



c coh {A, ^ r (A))a^ x0 (l + h(a, Q min )cosa) 



(31) 



and an analogous one for the spin contribution, h is the modulation ampli- 
tude and a is the phase of the Earth, which is zero around June 2nd, see 
Fig 3. The modulation amplitude would be an excellent signal in discrim- 
inating against background, but unfortunately it is very small, less than 
two per cent. Furthermore for intermediate and heavy nuclei, it can even 
change sign for sufficiently heavy LSP (see Fig. 4). So in our opinion a 
better signature is provided by directional experiments, which measure the 
direction of the recoiling nucleus. 



11 




Figure 3. The motion of the earth around the sun. The sun is moving with a velocity 
in the plane of the galaxy, which is perpendicular to the page. 




4.2. DIRECTIONAL RATES. 

Since the sun is moving around the galaxy in a directional experiment, i.e. 
one in which the direction of the recoiling nucleus is observed, one expects 
a strong correlation of the event rate with the motion of the sun. In fact 
the directional rate can be written as: 



Rdir — TT K 
Z7T 



c^A, fjL r (A))a^ x o(l + h m cos(a - a m tt)) 



(32) 
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TABLE 3. The parameters t, h, k, h m and a m for Qmin = 0. 
The results shown are for the light systems. +x is radi- 
ally out of the galaxy , +z is in the the sun's motion and 
+y vertical to the plane of the galaxy so that ((x,y,x) is 
right-handed. a m = 0, 1/2, 1, 3/2 implies that the maximum 
occurs on June, September, December and March 2nd re- 
spectively. 



type 


t 


h dir 


K 


h m 


a m 






+z 


0.0068 


0.227 


1 


dir 




+ (")x 


0.080 


0.272 


3/2(1/2) 






+(-)y 


0.080 


0.210 


0(1) 






-z 


0.395 


0.060 






all 1.00 

all 0.02 



and an analogous one for the spin contribution. The modulation now is h m , 
with a shift a m ir in the phase of the Earth a, depending on the direction 
of observation. k/(2tt) is the reduction factor of the unmodulated direc- 
tional rate relative to the non-directional one. The parameters k , h m , a m 
strongly depend on the direction of observation. We prefer to use the pa- 




0.5 1 1.5 2 2.5 3 07b 1 1^5 2 2^5 3 



Figure 5. The expected modulation amplitude h m for A — 127 in a direction outward 
from the galaxy on the left and perpendicular to the galaxy on the right as a function 
of the polar angle measured from the sun's velocity. For angles than 7r/2 it is irrelevant 
since the event rates are tiny. 

rameters k and h m , since, being ratios, are expected to be less dependent 
on the parameters of the theory. In the case of A = 127 we exhibit the the 
angular dependence of h m for an LSP mass of m x = lOOGey in Fig. 5. We 
also exhibit the parameters t, h, n, h m and a m for the target A = 19 in 
Table 3 (for the other light systems the results are almost identical). 
The asymmetry in the direction of the sun's motion [21] is quite large, 
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« 0.97, while in the perpendicular plane the asymmetry equals the modu- 
lation. 

For a heavier nucleus the situation is a bit complicated. Now the param- 
eters k and h m depend on the LSP mass as well, (see Figs 6 and 7). The 
asymmetry and the shift in the phase of the Earth are similar to those of 
the A = 19 system. 



0.07 
0.065 



0.055 




Figure 6. The parameter k as a function of the LSP mass in the case of the A = 127 
system, for Q m in = expected in a plane perpendicular to the sun's velocity on the left 
and opposite to the sun's velocity on the right. 



0.19 
0.18 




Figure 7. The modulation amplitude h m in a plane perpendicular to the sun's velocity 
on the left and opposite to the suns velocity on the right. Otherwise the notation is the 
same as in Fig 6. 



4.3. TRANSITIONS TO EXCITED STATES 

Incorporating the relevant kinematics and integrating the differential event 
rate dR/du from u m i n to u max we obtain the total rate as follows: 

fU m ax A f? 7/2 rUmax AT? 

R exc = / =(1 - ^f)du ,R gs = ^du (33) 

where u exc = J^^q and E x is the excitation energy of the final nucleus, 

Umax = (y/a) 2 - (E x /Q ) , y = v/upsilon and u min = Q min /Qo, Qmin 
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itilsp 



(GeV) 



Figure 8. The ratio of the rate to the excited state divided by that of the ground state 
as a function of the LSP mass (in GeV) for 127 1. We assumed that the static spin matrix 
element of the transition from the ground to the excited state is a factor of 1.9 larger 
than that involving the ground state, but the functions -Fii(it) are the same. On the left 
we show the results for Qmin = and on the right for Qmin = 10 KeV. 



(imposed by the detector energy cutoff) and u max = (y esc /a) 2 is imposed 
by the escape velocity (y esc = 2.84). 

For our purposes it is adequate to estimate the ratio of the rate to the 
excited state divided by that to the ground state (branching ratio) as a 
function of the LSP mass. This can be cast in the form: 

BRR ^* exc (ty ^ exci^exci Uumax)\\ ~\~ h exc {u exCl U max ) COS Ct] (34) 

S gs (0) *j S («mm)[l + h(u min ) cosa] 

in an obvious notation [21]). S gs (0) and S exc (0) are the static spin matrix 
elements As we have seen their ratio is essentially independent of super- 
symmetry, if the isoscalar contribution is neglected. For 127 I it was found 
[11] to be be about 2. The functions ^ are given as follows : 

Ay/a) 2 a („.) 

/ ir\^ F i^ u )^ a ^)-^yesc)]du (35) 



2 



/■Umax Q oi) 

y exc (u exc ,u max ) = / ^±F% c (u)(l-^- xc 



2 



[il)(a-\/u(l + u exc /u)) - ip{yesc)] du (36) 

The functions ip arise from the convolution with LSP velocity distribution. 
The obtained results are shown in Fig. 8. 



4.4. DETECTING RECOILING ELECTRONS FOLLOWING THE 
LSP-NUCLEUS COLLISION. 

During the LSP-nucleus collision we can have the ionization of the atom. 
One thus may exploit this signature of the reaction and try to detect these 
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m LS p -> (GeV) 



Figure 9. The the same as in Fig. 4 for the modulation amplitude h exc for the transition 
to the excited state. 



electrons [23]. These electrons are expected to be of low energy and one 
thus may have a better chance of observing them in a gaseous TPC detec- 
tor. In order to avoid complications arising from atomic physics we have 
chosen as a target 20 Ne. Furthermore, to avoid complications regarding the 
allowed SUSY parameter space, we will present our results normalized to 
the standard neutralino nucleus cross section. The thus obtained branch- 
ing ratios are independent of all parameters of supersymmetry except the 
neutralino mass. The numerical results given here apply in the case of the 
coherent mode. If, however, we limit ourselves to the ratios of the relevant 
cross sections, we do not expect substantial changes in the case of the spin 
induced process. 

The ratio of our differential (with respect of the electron energy) cross 
section divided by the total cross section of the standard neutralino-nucleus 
elastic scattering, nuclear recoil experiments (nrec), takes [23] the form: 

— = jE^I<M2m e T)| 2 

Puree 4 n£ 

/i^/i^^l^t + AMP^ (37) 



/ 1 2«[F(2 Mr «f)] 2 



where 



Pv -k J Pi + k 2 - 2kp x £ >1 

K = 5* , K = V-^ , a = p x .k, £ = q.K, 

Px Px 



£l = y^[l + 7^(^-7^ — 1)] an d 2^p x £ = 2fi r v( is the momentum q 
transferred to the nucleus and the momentum q transferred to the nucleus 
and F(q) is the nuclear form factor. The outgoing electron energy lies in the 
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Figure 10. Shown on the left is the differential rate, divided by the total rate associated 
with the nuclear recoils, as a function of the electron energy T (in keV). Each atomic 
orbit involved in the target 20 Ne is included separately. The full line, the short-dashed 
line and the long-dashed line correspond to the orbits Is , 2s and 2p respectively. Shown 
on the right is the ratio of the total rate for the novel process divided by that of the 
standard process as a function of the electron threshold energy, assuming zero threshold 
energy for the standard process. This ratio may increase if such a threshold is included. 



range < T < ^-T x — e n £. Since the momentum of the outgoing electron is 
much smaller than the momentum of the oncoming neutralino, i.e. K ~ 1, 
the integration over £i can be trivially performed. 

We remind the reader that the LSP- nucleus cross-section a nrec takes 
the form: 

Onrec = (-^-) 2 ^V p / 2<%[F(2fl r v£)] 2 (38) 

In the case of 20 Ne he binding energies and the occupation probabilities are 
given by [23]: 

e ne = (-0.870, -0.048 , -0.021) , p nt = (2/10, 2/10, 6/10). (39) 

in the obvious order: (ls,2s,2p). In Fig. 10 we show the differential rate of 
our process, divided by the total nuclear recoil event rate, for each orbit as 
well as the total rate in our process divided by that of the standard rate as 
a function of the electron threshold energy with threshold energy in the 
standard process. We obtained our results using appropriate form factor 
[13]. 

From these plots we see that, even though the differential rate peaks 
at low energies, there remains substantial strength above the electron en- 
ergy of 0.2 keV, which is the threshold energy for detecting electrons in a 
Micromegas detector, like the one recently [24] proposed. 5. Conclu- 
sions 

Since the expected event rates for direct neutralino detection are very 
low [6, 9], in the present work we looked for characteristic experimental 
signatures for background reduction, such as: 

Standard recoil experiments. Here the relevant parameters are t and 
h. For light targets they are essentially independent of the LSP mass [21], 
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essentially the same for both the coherent and the spin modes. The modu- 
lation is small, h « 0.2%, but it may increase as Q m in increases. Unfortu- 
nately, for heavy targets even the sign of h is uncertain for Q m i n = 0. The 
situation improves as Q m in increases, but at the expense of the number of 
counts. 

Directional experiments [19]. Here we find a correlation of the rates 
with the velocity of the sun as well as that of the Earth. One encounters re- 
duction factors n/2ir, which depend on the angle of observation. The most 
favorable factor is small, « 1/4-7T and occurs when the nucleus is recoiling 
opposite to the direction of motion of the sun. As a bonus one gets mod- 
ulation, which is three times larger, h m =« 0.06. In a plane perpendicular 
to the sun's direction of motion the reduction factor is close to 1/127T, but 
now the modulation can be quite high, h rn « 0.3, and exhibits very in- 
teresting time dependent pattern (see Table 3. Further interesting features 
may appear in the case of non standard velocity distributions [18]. 
Transitions to Excited states. We find that branching ratios for tran- 
sitions to the first excited state of 127 1 is relatively high, about 10%. The 
modulation in this case is much larger h exc ~ 0.6. We hope that such a 
branching ratio will encourage future experiments to search for character- 
istic 7 rays rather than recoils. 

Detection of ionization electrons produced in the LSP collision 

Our results indicate that one can be optimistic about using the emitted 
electrons in the neutralino nucleus collisions for the direct detection of the 
LSP. This novel process may be exploited by the planned TPC low energy 
electron detectors. By achieving low energy thresholds of about 0.25 keV, 
the branching ratios are approximately 10 percent. They can be even larger, 
if one includes low energy cutoffs imposed by the detectors in the standard 
experiments, not included in the above estimate. 

As we have seen the background problems associated with the proposed 
mechanism are not worse than those entering the standard experiments. 
In any case coincidence experiments with x-rays, produced following the 
de-excitation of the residual atom, may help reduce the background events 
to extremely low levels. 
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